Implementacja w MATLABie schematu alpha dla problemu transportu ciepła
Poniżej przedstawiam kod MATLABa wykonujący symulacje dla schematu alfa dla dwuwymiarowego problemu transportu ciepła.
Wykonanie kodu możliwe jest również w darmowym środowisku Octave.
Pobierz kod lub zob. Załącznik 7.
W linii 621 podajemy rozmiar kroku czasowego \( dt=0.01 \), w linii 622 podajemy parametr
\( alpha=0.5 \), pamiętając że alpha=0 to metoda explicite Euler, alpha=1 to metoda implicite Euler, alpha=0.5 to metoda Crancka-Nicolson.
W linii 623 podajemy liczbę kroków czasowych
\( K = 20; \)
Kod można uruchomić w darmowym środowisku Octave.
Kod uuchamia się otwierając go w Octave oraz wpisując komendę
\( heat\_time \)
Podczas działania kod wypisuje kolejne kroki czasowe
Iter 1, t = 0.010000
Iter 2, t = 0.020000
Iter 3, t = 0.030000
...
W każdej chwili czasowej w katalogu bieżącym kod generuje plik out_*.png, np.
out_0.png
out_1.png
out_2.png
...
zawierające rozwiązania z poszczególnych kroków czasowych.